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We present a scheme of positioning a single electron spin with sub-nanometer resolution through 
multiple nitrogen-vacancy centers in diamond. With unwanted noise suppressed by dynamical de¬ 
coupling, the spin coherence of each center develops characteristic oscillations due to a single electron 
spin located 20 nm away from the centers. We can extract the position information from the 
characteristic electron spin-coherence oscillations of each center. This scheme is useful for high- 
resolution nanoscale magnetometry. 

PACS numbers: 03.65.Yz, 76.30.Mi, 76.60.Lz 


I. INTRODUCTION 

Magnetic resonance spectroscopy of ensembles of spins 
has wide applications in many fields such as material sci¬ 
ence, analytical chemistry, structural biology [l|. Push¬ 
ing the sensitivity of conventional magnetic resonance 
imaging to the single spin level would enable even more 
important application in single molecule structure anal¬ 
ysis ii- However, it is extremely difficult to detect the 
weak magnetic field from a single electron or nuclear spin 
by conventional magnetic resonance spectroscopy [J, l5| . 
The magnetic resonance force spectroscopy has been de¬ 
veloped to improve the sensitivity to a few spins but with 
requirements of low temperature and high vacuum [^, Q , 
which limit its application to biological systems. 

Recently single spin detection in ambient conditions 
has been made possible by employing the atomic-scale 
defect in diamond—the nitrogen-vacancy (NV) center [1- 
El —as a magnetic sensor. The NV center has a para¬ 
magnetic spin-triplet ground state that can be optically 
initialized and read out [l^, [l^ and has a long spin coher¬ 
ence time up to milliseconds even at room temperature 
M- With recent development in diamond nanofabri¬ 
cation [^, the NV center in high-purity diamond has 
been used as a magnetic sensor to detect single nuclear 
spins (l6l - [T8l| and nuclear spin clusters [l^ in diamond. 
The shallow NV centers below the diamond surface have 
also been used to sense single electron spins in single 
molecules [1^ [2l| and electron spin labels in single pro¬ 
teins on the surface. The high-resolution positioning 
of single spins has potential applications in revealing the 
structure and dynamics of single molecules . 

A widely-used sensing scheme is the quantum sensing 
based on dynamical decoupling (DD) [23,[il| based quan- 
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turn sensing. The key idea is to identify the characteris¬ 
tic oscillations caused by the target spin imprinted onto 
electron spin coherence of the NV center when the noise 
from the target spin is resonantly amplihed by DD con¬ 
trol of the NV electron spin [2^ . The DD control can 

also suppress the background noise and prolong the co¬ 
herence time of NV electron spin. However, the position 
of the target spin cannot be uniquely determined by one 
experimental setting [illlS- Instead, the measurement 
has to be repeated for various magnetic field orientations 
or for different scanning positions of the sensor to fully 
position the target spin. Such a requirement restricts the 
speed of scanning and limits the potential applications in 
single-molecule structure imaging. 

In this paper, we propose to use multiple NV centers 
in a diamond probe tip for positioning single spins with 
sub-nanometer resolution [Fig. [l][a)]. This scheme is 
analogous to the multi-satellite positioning systems. For 
a given NV center spin, the coherent oscillations provide 
partial information on the distance and direction of the 
target single spin. The full position information of the 
target spin can be obtained by the oscillations of the 
multiple NV center spins. This scheme would greatly 
improve the speed of detection compared with the scheme 
using one NV center [TMllil , since there is no need 
of varying the magnetic field or nanometer-step scanning 
the probe [T^ . 


II. THEORETICAL MODEL 

As shown in Fig. [TJa), the proposed setup to detect 
a remote single electron spin contains three NV centers 
in a diamond dip. The spin Hamiltonian of the whole 
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FIG. 1: Fingerprint features in NV center spin coherence due to coupling to a single electron spin, (a) Schematic illustration of 
the single electron spin detection based on three NV centers embedded in a diamond tip. The NV axes of the three NV centers 
are all chosen along the crystal axis [111]. The close-up shows the atomic structure of the NV center, (b) The relative position 
between a single NV sensor spin and the target electron spin is denoted by the distance R and the zenith angle 6 relative to 
the [111] axis, (c) NV center spin coherence under CPMG-30 control as a function of time and the zenith angle of the target 
spin at a distance of i? = 5 nm. (d) NV center spin coherence under CPMG-30 control as a function of time and the distance 
R of the target spin with a zenith angle 0 = 30°. The strain-induced transverse anisotropy parameters are ea = 3 MHz, eb = 2 
MHz, £c = 4 MHz and a magnetic held B — 0.1 Gauss is applied along the NV axis. The target electron spin is assumed to 
have the same gyromagnetic ratio as the NV electron spins (7e = 7nv). The calculations in (c) and (d) are performed for one 
of the three NV centers (NV-A) with only the coupling between NV-A and the target spin considered. 


system is 

H =HtaY + ^ + ^^bath + ^ ^^NVi-NVj 

i i,j 

—tar E — bath; (1) 

i i 

where the first three terms denote the Hamiltonians of 
the isolated target electron, the NV center sensors and 
the nuclear spin bath correspondingly, and the last three 
terms represent the interaction between different NV cen¬ 
ter sensors, the coupling of the NV center sensors to the 
target electron spin, and the coupling of the NV center 
sensors to the nuclear spin bath in diamond. We assume 
that the Hamiltonian of the spin-1/2 target electron spin 
is iTtar = ~ 7 eB • Sg with 7 e and Se denoting the gyro- 
magnetic ratio and spin operator of the target spin re¬ 
spectively, and the Hamiltonian of the Hh NV center is 


where Si is electron spin operator of the zth NV center, 
A = 2.87 GHz is the zero-field splitting of the NV centers, 
Bi denotes the local strain-induced transverse anisotropy 
for the ith NV center, 7nv=-2.8 MHz Gauss“^ is the gy¬ 
romagnetic ratio of NV electron spin, B is the magnetic 
field, and the z axis is the intrinsic electron spin quantiza¬ 
tion axis of the ith NV center (pointing from the nitrogen 
to the vacancy). In this paper, 

With a rather strong magnetic field applied along 
2 ; axis, the local strain terms in Eq. m can be ne¬ 
glected and the ith NV center has three spin eigenstates 
{|0)i, |±l)i}. However, in the weak magnetic field regime, 
the strain-induced transverse anisotropy breaks the de¬ 
generacy of the I ± l)i states in zero magnetic field (the 
eigenstate |0)i is almost unperturbed since A,. ^ £i) and 
results in renormalized eigenstates |±)i (see Ap¬ 

pendix A). In the new eigenstate basis {|0)i, |-|-)i, |—)i} 
of the ith NV center, the full Hamiltonian in Eq. (1) can 
be recast into the pure-dephasing form as (see Appendix 

B) 

H, = hI+^\+U+\ + Ht^\-U-\ + |0).(0|, (3) 


iJNv. = HSt? + eAiSt? - - 7nvB • S„ (2) 
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=HtaT + ^2 -NVfc + ^^bath 

j,k^i 

+ ^ (^^NVj + H^Vj -t 
3^'i- 


-bath } 


(5) 


where wnv = |7nv|-S is the Larmor frequency of the NV 
electron spin, hi = hf + + hf is the noise field for the 

ith NV center, hf = {z- A^^e) ■ Sg is the dipolar magnetic 
field produced by the target electron spin, ' 

'Im is the nuclear Overhauster field with denot¬ 
ing the nuclear spin operator of the mth nuclear spin in 
the bath, hf = ' ^i,j)' is the dipolar magnetic 

field produced by the other NV centers, Ai e,Ai ,„,Ai j 
are the dipolar interaction tensors for the coupling of the 
ith NV electron spin to the target spin, the mth nuclear 
spin and the jth NV electron spin correspondingly, and 
-f^bath = 7nB • Im “f where 'Jn 

is the gyromangetic ratio of nuclear spins and 'Dmn is the 
nuclear-nuclear dipolar interaction tensor. 

As shown in Eq. o, the renormalized eigenstates 
|±)i have eigenenergies ±\/ef + Wnv renormal¬ 

ized eigenstates |±)i, the coupling strength of the ith 
NV center to the target spin, the other NV centers and 
the nu clear spin bath are all reduced to a factor of 
WNv/\/£i + Moreover, The last term in Eq. (4) 

mainly induces long-range interaction within the nuclear 
spin bath, but its effect on the sensor coherence can be 
largely suppressed for not too small magnetic field and 
large strain parameters. In this paper we consider the 
single transition (|0) O |-l-)i) of all the NV sensors. 

The NV center spin decoherence in high-purity dia¬ 
mond is mainly caused by the hyperfine interaction with 
the nuclear spins [3 lUl- At finite temperature, 
the random orientations of the nuclear spins result in 
local field fluctuation (thermal fluctuation) that can be 
eliminated by spin echo. The internal dynamics of the 
spin bath also induces the dynamical quantum fluctua¬ 
tion [s^l , which cannot be completely removed by dynam¬ 
ical decoupling. To detect a single electron spin by the 
NV center, we should suppress the thermal and quantum 
noise caused by the nuclear spin bath while amplifying 
the effect of the target electron on the NV center spin co¬ 
herence. The Carr-Purcell-Meiboom-Gill (CPMG) con¬ 
trol [13, m, 1131 meets both requirements, since it can 
simultaneously increase the coherence time of the sensor 
spin while selectively amplifying the noise at a specific 
frequency. For the V-pulse GPMG (CPMG-V) sequence, 
the NV electron spin is flipped at time tk = (2fc— 1 )/ (27V) 
with t being the total evolution time and k = 1,2, ■ ■ ■ ,7V. 
The target spin can be detected if the NV spin coherence 
profile develops characteristic oscillations caused by the 
target spin. 


The oscillation features in the NV sensor spin coher¬ 
ence caused by the target spin can be understood as fol¬ 
lows. In the classical noise picture, the precession of the 
single electron spin near the NV sensor gives rise to addi¬ 
tional peak structures 12711 i n the smooth noise spectrum 
of the nuclear spin bath [3Ij , corresponding to a series of 
coherence dips in the NV sensor spin decoherence pro¬ 
file. When the magnetic field is much stronger than the 
dipolar interaction with the target spin, the peak in the 
noise spectrum is approximately located at the electron 
spin Larmor frequency uje = \'Je\B. For GPMG-TV con¬ 
trol, the coherence dips of the 7th NV sensor spin caused 
by a remote target electron spin occur approximately at 
times 


7r(2g- 1)V 
We -f XiAf^^/2 


( 6 ) 


where q = 1, 2, • • • denotes the dip order. In the following 
we will always consider the first-order coherence dip {q = 
1). The sensor coherence dip depth as a function of the 
GPMG pulse number can be derived in the weak sensor- 
target coupling regime as 


Adip(A^)«cos 


\We -I- AiAf g/2y ’ 


(7) 


where Xi = wnv/ is the renormalization fac¬ 

tor for the 7th NV sensor, Afg, A^g are the com¬ 
ponents of the dipolar interaction with the target 
spin parallel and orthogonal to the 2 axis respec¬ 
tively. For a small pulse number, the above equa¬ 
tion agrees with that from the semiclassical noise model 
7-dip(A^)«exp[-V2(A,A^g)72(a;e -f A,Af_g/2)2] (see Ap¬ 
pendix C and D for details). 


III. DETECTION SCHEME 

To detect the target spin with high spatial resolution, 
the target spins in different positions should have distin¬ 
guishable fingerprint oscillations imprinted onto the sen¬ 
sor spin coherence. This implies that the applied mag¬ 
netic field should be comparable to the dipolar interac¬ 
tion between the sensor spin and the target spin. For 
detecting a single electron spin 5~I0 nm away, the mag¬ 
netic field should be about 0.2~0.05 Gauss. The renor¬ 
malization factor is about 0.18~0.05 (the strain-induced 
transverse anisotropy ^ 3 MHz), which means the 
effective dipolar interaction is about ten times smaller 
than the magnetic field. To increase the detection res¬ 
olution, we should use the GPMG control with a large 
number of pulses. 

The renormalized dipolar interaction between the 
7th sensor spin and the target spin has components 
[A|j, A7 = [1 — 3cos^0, 3sin0cos0] (where /xq 

is the vacuum permeability) that are determined by the 
distance R and the zenith angle 6 between the displace¬ 
ment and the [III] axis [Fig. (Hb)]. Figure [T](c) and[T](d) 
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FIG. 2: Oscillations caused by the target spin on the NV cen¬ 
ter spin decoherence with the nuclear spin noise suppressed 
by dynamical decoupling, (a) Decoherence of NV-A center (in 
the absence of the target spin and other NV centers) caused 
by the *^^0 spin bath with a natural abundance (1.1%) un¬ 
der CPMG-30 control, (b) Spin coherence of three NV cen¬ 
ters (with the presence of the target spin) under GPMG-30 
control (circles) which matches the contribution solely from 
the target electron (solid lines) located at {Ra = 7.46 nm, 
Oa = 19.56°,i?s = 8.72 nm, Ob = 33.92° and Rc = 8.83 nm, 
6c = 35.03°). The parameters are the same as those in Fig. 

m 


show the different coherence oscillations, from which we 
can see that an atomic-scale change of the position of the 
target spin would have noticeable influence on the time- 
domain coherence features of the NV sensor spin. It is 
worth noting that the contour plot in Fig. [ijc) is mirror 
symmetric about the 6 = 90° line owing to the symmetry 
of the dipolar interaction. We may use the coherence fea¬ 
tures to determine the distance and zenith angle of the 
target spin. 

To fully determine the three-dimensional position of 
the target spin, we propose to use multiple NV centers as 
the sensor spins, in analogue to the multi-satellite global 
positioning systems. We note that recently strongly cou¬ 
pled NV pairs below the diamond suface with the cou¬ 
pling strength up to 53 kHz (corresponding to the dis¬ 
tance about 10 nm) and an average depth of 15 nm have 
been successfully generated by implanting ionized nitro¬ 
gen molecules (^®N^) into diamond [d^. In principle, 
multiple closely spaced NV centers can be generated by 
implanting N 3 and N 4 molecules into diamond 

We consider three NV sensor spins (A, B, C) embed¬ 
ded near the surface of a diamond tip. The three NV 
sensors can be separately addressed by their different en¬ 
ergy splittings between |0)i and |-|-)i due to the difference 
in the local strain that usually occurs in real diamond 
[^ . In the simulation, we assume the strain-induced 
transverse anisotropy parameters of the three NV sen¬ 
sors to be £a = 3 MHz, sb = 2 MHz and ec = 4 
MHz. The differences in the local strain also help sup¬ 
press the exchange interaction between the NV sensor 
spins in the multi-spin sensor, which simplifies the data 
analysis for the positioning. We may use the CPMG 
control with a specific frequency to induce the spin tran¬ 
sition in one NV sensor while leaving the states of the 
other NV sensors unchanged. We get three sets of co¬ 


ordinates {Ra, Sa),{Rb,Sb), {Rc,dc) by comparing the 
coherent oscillations due to the target spin with the dip 
features in the spin decoherence of the three NV sensors. 
Therefore the three-dimensional coordinates of the target 
spin are determined. 

For simplicity, we suppose the three NV sensors are all 
along the crystal axis [111]. If they have different direc¬ 
tions in the diamond lattice, the three sets of coordinates 
obtained by the detection scheme are with respect to the 
three different local coordinates of the three sensors, and 
need to be transformed into the same coordinate for po¬ 
sitioning the target spin, and the effective Zeeman split¬ 
tings relevant to the renormalization factors Xi of the 
NV sensors (projections of the magnetic field to the NV 
sensors’ axes) are different for different sensors. 

The NV sensors are separated from each other by 6 7 

nm with the dipolar interaction between them less than 
0.1 MHz, while the energy cost of pairwise flip-flop pro¬ 
cesses (e.g. |0)i|-|-)j ■<->■ |+)i|0)j) between two NV sen¬ 
sors is approximated as |ei — e^j ^ MHz, which is much 
stronger than the dipolar interaction between them. As 
a result, the flip-flop processes between different NV elec¬ 
tron spins are highly suppressed [1^. Moreover, the sin¬ 
gle spin dynamics of the other NV sensors on the deco¬ 
herence of a specific NV sensor is very small due to the 
much larger noise frequency from the ith NV sensor (cor¬ 
responding to the transition |+)i_0 j— )i) compared with 
target spin Larmor frequency [30|. Therefore for a spe¬ 
cific NV sensor, the effects of other sensors on the sensor 
decoherence can be neglected. 

We solve the pure dephasing problem with inclusion 
of the interactions between the NV sensors, the coupling 
to the target electron spin, and the coupling to the nu¬ 
clear spin bath, by adopting the well-established cluster- 
correlation expansion method (CCE)[33,[3l[3lliIl- The 
key idea is that the qubit coherence can be expressed as 
the product of cluster correlations. In real calculations, 
it often suffices to truncate the expansion up to the min¬ 
imum size M of the clusters (CCE-M) to get converged 
results. It is clear from Fig. Hla) that the spin coherence 
of a single NV sensor in nuclear spin bath with a 
natural abundance 1 . 1 % can be well protected for t < 1 
ms under CPMG-30 control. 

Apart from the nuclear spin bath, the electron 
spin bath on the diamond surface cause additional de¬ 
coherence for the near-surface NV centers In 

recent experiments, stable and well-behaved NV centers 
located with depth ranging from 1 to 10 nm from the 
diamond surface have been observed and the coherence 
time T 2 longer than 100 p,s has been achieved for the 
depth of 5 nm from the diamond surface . The 

coherence time of the near-surface NV electron spin can 
be further increased by high temperature annealing [d^ 
and etching away the diamond surface [d^. As shown 
in Figs. [T][c) and[T][d), the typical timescale of the fin¬ 
gerprint oscillation of a target spin is much shorter than 
the NV electron spin coherence time, so the nuclear spin 
bath would have negligible effect on the positioning. The 
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FIG. 3: Numerical simulation for positioning a target spin. 
The red lines represent the spin coherence obtained from the 
detection while the black solid lines represent the matched 
oscillation patterns in the fingerprint library (the discrete 
steps are taken as dR = 0.02 nm, d6 = 0.2° here). The 
position and depth of the first dip in the spin coherence os¬ 
cillations are used as the criteria for matching. The esti¬ 
mated locations of the electron relative to the three NV cen¬ 
ters are (a) Ra € [7.40 ~ 7.50 nm], 9a G [18.4° - 20.0°], 
(b) Rb G [8.64 - 8.76 nm], Ob G [32.2° ~ 34.2°[ and (c) 
Rc G [8.64 ~ 8.96 nm], Be G [33.0° ~ 36.0°]. The ex¬ 
act position of the target spin is such that [Ra = 7.46 nm, 
6A = 19.56°, i?fl = 8.72 nm, Ob = 33.92° and Rc = 8.83 
nm, 6c ~ 35.03°). The green lines represent the oscillation 
patterns if the target spin is moved about 0.6 nm away from 
its original position. The parameters are the same as those in 

Fig.m. 

other NV sensors also have negligible influence on the 
sensor spin coherence due to the different strain factors 
of different NV sensors, therefore the spin coherence of 
the NV sensors matches the contribution solely from the 
target electron spin [Fig. EKb)]. 

Figure [3] shows the numerical simulation for position¬ 
ing a target spin. First, we establish a fingerprint library 
to store the positions and depths of the dips or peaks in 
the decoherence patterns caused by a single electron spin 
located in a large range (5 nm<i? < 30 nm, 0 < 0 < 7r/2). 
R and 0 are discretized with resolution (ii?=0.02 nm and 
d0=O.2°. Second, we put the diamond tip containing 
three NV centers close to the target spin and get the re¬ 
spective spin coherence profiles of the three NV centers 
under CPMG-30 control. Third, by matching the sensor 
spin coherence oscillations to the fingerprint library, the 
position of the target spin is restricted to a small spatial 
region intersected by the small ranges of the three sets 
of parameters (Ri,0i). The smaller the volume of the 
intersected region, the higher the resolution in detecting 




FIG. 4: Positioning of a single spin about 10~20 nm away, 
(a) NV center spin decoherence caused solely by the target 
spin under GPMG-100 control as a function of time and the 
distance R with a zenith angle 9 — 30°. (b) Spin coherence 
of three NV centers under CPMG-100 control (dashed lines 
with circles) which matches the contribution solely from a 
single electron spin (solid lines) located at (Ra = 13.13 nm, 
6a = 23.60°, Rb = 12.57 nm, 9b = 13.35°, = 14.73 nm, 

9c = 33.88°). Here the additional narrow dips in the detected 
signal of a specific NV sensor (dashed lines with circles) is 
caused by the other NV sensors. The parameters are the 
same as those in Fig. [T] 

the target spin. In the simulation, the resolution of less 
than 0.3 nm is achieved. 

Now we discuss how to extend the detection range of a 
target spin. The strain-induced transverse anisotropy of 
the NV sensor leads to the energy splitting of the | ± l)i 
states, and the sensor in the renormalizd eigenstates |±)i 
have a smaller effective magnetic moment than that of 
the sensor in states |±l)i. Consequently, the decoherence 
caused by the target spin would be suppressed, which is 
disadvantageous to the positioning. To increase the de¬ 
tection range, the strain-induced transverse anisotropy 
should be decreased. Moreover, since the dipolar inter¬ 
action between the NV sensor spin and the target spin 
decreases rapidly as the distance between them increases, 
the magnetic field should be decreased correspondingly 
to be comparable to the interaction between the sensor 
spin and target spin. However, the renormalized dipolar 
interaction decreases almost linearly with the magnetic 
field. Therefore, the detection range cannot be increased 
by decreasing the magnetic field. In this case, we can in¬ 
crease the number of CPMG pulses to further amplify the 
effect of slight difference of the dipolar interaction on the 
time and depth of the coherence dip. As shown in Fig. 
m under CPMG-IOO control, the detection range of the 
multi-spin sensor can be extended to ^20 nm. At such a 
large distance, we can still achieve sub-nanometer reso¬ 
lution in the simulation. Moreover, with the extension of 
the detection range, we can choose the near-surface NV 
centers with larger depth (> 10 nm) so that the surface 
decoherence effects can be greatly reduced (42l - l4^ . 

To increase the CPMG pulse number with a fixed mag¬ 
netic field requires a longer coherent evolution time for 
both the sensor and target spin. For the case in Fig. IH 
the target spin should have long coherence time greater 
than 100 ^s. However, this requirement can be relaxed 
by increasing the external magnetic field according to Eq. 
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(| 6 |) . A key requirement to realize nanoscale resolution in 
single spin detection in our scheme is that the external 
magnetic field should be comparable to the dipolar inter¬ 
action between the sensor and target spin. If the hyper- 
fine interaction is increased by lowering the strain factors, 
which may be realized by choosing deeper NV centers 
or tuning of the strain factors by electric fields , then 
the magnetic field can be increased without decreasing 
the spatial resolution. For example, if the strain factor 
Ei of the fth NV sensor is decreased from 3 MHz to 0.3 
MHz, the magnetic field can be increased from 0.1 G to 1 
G with the ratio of We to /V^^i + unchanged, 

so the minimal coherent evolution time of the target spin 
is decreased to about 5 ^s for CPMG-30 and 18 for 
CPMG-100. 

In Ref. [^ . a scanning magnetometer containing one 
NV center was employed to realize nanoscale magnetic 
imaging of another target electron spin located about 50 
nm away. There the sensor and target spins were manip¬ 
ulated synchronously by periodic DD control (up to 100 
pulses), to constructively accumulate the phase shifts in¬ 
duced on the sensor by the target spin. In our scheme, 
only the sensor spin is manipulated by DD control, and 
the sensor spin coherence dip is caused by the intrinsic 
quantum dynamics of the target spin. 


IV. CONCLUSIONS 

In summary, we have proposed and numerically 
demonstrated that atomic-scale positioning of single elec¬ 
tron spins can be achieved by using multiple NV centers 
as the sensors. Each NV sensor spin works independently 
and gathers information about the spatial range of the 
target spin. By integrating the information provided by 
multiple NV sensor spins, the position of the target can 
be accurately determined. The scheme, without requir¬ 
ing spatial scanning or varying the magnetic field direc¬ 
tion, may provide an approach to fast and scalable posi¬ 
tioning of single spins with sub-nanometer resolution. 
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Appendix A: Diagonalization of the NV center 
Hamiltonian 


Finally, we discuss some details of the possible experi¬ 
mental realization of our detection scheme. To accurately 
position a target spin, we have to accurately determine 
the relative locations, strain parameters and orientations 
of the multiple NV sensor spins. Recently the depth of 
shallow NV centers in diamond can be determined with 
^ 1 nm resolution by detecting the nuclear magnetic res¬ 
onance signal of a proton nuclear spin bath placed on the 
surface [51| . The strain parameters of the NV sensors has 
also been determined with ~ 2 kHz accuracy by measur¬ 
ing ODMR spectra of NV sensors at zero or very weak 
magnetic field ■ The orientation of a NV sensor can be 
determined by measuring the optically detected magnetic 
resonance spectra (ODMR) in a relatively large magnetic 
field , since the sensor electron spin Zeeman splitting 
would be different for different orientations. Although we 
assume that all the three NV sensors have the same ori¬ 
entation, the detection scheme can still work for multiple 
NV sensors of different orientations by just transforming 
between different local coordinates of the different NV 
sensors when performing data analysis of different sensor 
signals. Moreover, the control errors of finite DD pulses 
can accumulate for a large DD pulse number and degrade 
the DD performance. In experiments, the XY 8 -fc pulse 
sequence rather than CPMG pulse sequence can largely 
suppress the pulse errors [^, . Recently the nuclear 

magnetic resonance of multiple nuclear species has been 
successfully detected by using a shallow NV sensor un¬ 
der up to XY8-40 control (equivalent to GPMG-320 in 
theory) [H^. 


The Hamiltonian of the ith NV center in Eq. © can 
be written in a matrix form in the basis of the eigenstates 
{| + 1)», |0)j, I - 1)J as 

/ A -1 - WnV 0 Ei \ 

Hnv, = 0 0 0 , (Al) 

y Ei 0 A — wnv / 

Diagonalization of this matrix gives the new basis of 
eigenstates {|+)i, |0)i, |—)i} with eigenenergies {A + 
0, A — -^/e^ -I- a;^y}. The new eigenstates 
are 


i±)i 


I 


^WNV ± 



+ l)z + £i 



where C± = y^ 2 -f w^( y'ef -f ± wn\0 • 


Appendix B: Derivation of the renormalized 
pnre-dephasing Hamiltonian 


Since the zero-field splitting of the fth NV sensor is 
much larger than its dipolar coupling to the target spin, 
the other NV sensors and the nuclear spin bath, the state 
| 0 )i can be safely assumed to be decoupled from the states 
I ± I)i. Then the effective Hamiltonian in the subspace 
{| -I- l)i, I - I)i} becomes 


H,=A + 


WnV + hi Ei 

£i —WnV “ hi 


H, 


(0) 


(Bl) 
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The diagonal form of this Hamiltonian is obtained as 


Hi — \Jef + (wnv + ^i)^(l+)i(+l “ l“)j(~l) + h\ ^(d32) 

where we have dropped the zero-field splitting A for sim¬ 
plicity. Since the target spin, NV sensors and nuclear 
spin bath are generally unpolarized at room tempera¬ 
ture, the mean value of the noise field hi is much smaller 
than y/sj + Hamiltonian can be expanded 

in the second order of hi as 

Hi = hI+'>\+)i{+\+hI-^\-)i{-1 (B3) 


with 


m 


(±) 





Now we can write the Hamiltonian in the full Hilbert 
space {|-l-)i, |0)i, |—)i} and get the pure-dephasing Hamil¬ 
tonian in Eq. 


Appendix C: General semiclassical noise model for 
qnantum sensing 


We consider a general Hamiltonian of a spin-1/2 sensor 
coupled to a target spin cluster [27|, 


H = S,P + Ho, 


(Cl) 


with 


d 

Ho=YiEn\n){n\, (C2) 

n—1 

/3=i^(^m„|m)(n|-bH.c.), (C3) 

m^n 


where P is the noise operator for the sensor, Hq is the free 
Hamiltonian for the target spin cluster, d is the number 
of eigenstates of the target spin cluster, Pmn = (w|/3|n). 
In the interaction picture set by Hq, the time-dependent 
noise operator is P{t) = ^ noise cor¬ 

relation function is just 

C{t) = mpm = i ^ (C4) 

m,n 


where = Em - En and (•••) = Tr[po • • • ] with po = 
^(^^1 |TO){m|. The noise spectrum is the Fourier 
transform of the noise correlation function , 


C{oj) = I \Pmn\'^S{u} - UJmn)- 

(C5) 


It is easy to demonstrate that C{t) = C{—t) and C{uj) = 
C{—uj). Typically DD control (consisting of a sequence 
of TT flips at times {ti, t 2 ■ • ■ tAr} for the sensor evolution 
from 0 to t) is applied to the sensor spin to selectively 
amplify a specific noise frequency. With the Gaussian 
noise approximation, the sensor spin decoherence under 
DD control is [3l| 


L{t) = exp 


= exp 



J dtidt2C{ti 




1 

2 





(C6) 


where fit) = (—1)^ for [tk,tk+i] is the DD modulation 
function (to=0, tAr+i = t), F{uj,t) = w|/J/(t')dt'p = 

is the DD filter function. 

For the A-pulse CPMG control, tk = (2fc — l)/2iV 
(fc = 1, 2, • • • , N) with the pulse delay 2r = t/N and the 
filter function is [3l| 


E{uj,t)= 


4sin2(^l 


4sin 


\4NJ 


cos (^f) cos i(^)|,oddiV,, . 

sin (f) cos-1 ^i^jl^g^gn jvrG 


If the CPMG pulse delay matches one noise frequency 
\LUmn\, that is 2rdip = (2g - l)7r/|wmn| or tdip = 7r(2g - 
l)N/\u:mn\ {q = 1, 2, • • •), the noise with this frequency 
is amplified by CPMG control and the sensor spin coher¬ 
ence shows sharp dips. Since the filter function at the 
coherence dips is F'(a;m„,tdip) = 2A, the sensor coher¬ 
ence dip depth corresponding to \u}mn\ as a function of 
CPMG pulse number is 


idip(^) 



4N^\Pmu? \ 

) 


(C8) 


The semiclassical noise model is valid for the weak sensor- 
target coupling regime (|/3mn| l^mnl) and relatively 
small DD pulse number. Moreover, the sensor coherence 
dip depth from the semiclassical noise model is always 
positive while the one from exact quantum model can be 
negative [T^ . 


Appendix D: Sensing a single target electron spin by 
the NV sensor spins 

The Hamiltonian for the ith NV sensor coupled to a 
remote target electron spin is 


H = {\ihl 
XM 


■CTj -l- 


Htax) |+)i(+| 

XiK 

- H- rifai 


i7tar|0),(0| 


(DI) 


where erf = |-l-)i(-l-| — |0)i(0| is the pauli operator in the 
subspace {|0)i, |-|-)i}. A* = WNv/y/ej + is the renor¬ 
malization factor, hf = (z • ■ Sg = e • Sg and 

Htai = ^eSg- In comparison with the general Hamil¬ 
tonian in Eq. the noise operator for the ith NV 


— OO 




















center is 


/3 = A,/i® = A,A,.e-Se, (D2) 

and the free Hamiltonian of the target spin is 

Ho = X,ht/2 + Htar = + u:,Sl (D3) 


Since ||/3|| <C ||-Ho||, we adopt the interaction picture, 
in which the noise operator is defined as 

/3(t) = e^Hotp^-iHot 

r 


5' + g»(‘^e*-“) _|_ 


(D8) 


1. Semiclassical noise model 

In the weak-coupling regime (||/3|| <C ||-Ho|| or |Ai_e| <g; 
We), the noise spectrum of the target electron spin can 
be simplihed as 

C(w) ^ - w;) + <5(w + w')], (D4) 


where w' and are the same as those defined in the last 
part, = Sg ± iS^, and a = arctan(A^g/Af g). Here 
we have made an approximation Hq ~ (A^Af g/2 + a;e)S'|. 
Then and can be written as 


where w' ss X^A\ J2 + and A^g = ^(Af_g)2 + (Af_g)2 

(the zero-frequency noise has no contribution to the sen¬ 
sor decoherence under DD control and is neglected here). 
We can see that the noise frequency of the target spin is 
the Zeeman frequency plus a shift caused by the hyper- 
fine interaction with the NV sensor spin. So the sensor 
coherence dip occurs at times tdip = 7J‘(2g — 1)A/Wg, and 
the coherence dip depth as a function of the CPMG pulse 
number is 


where /(t) = (—1)^ for \tk,tk+i\ is the DD modulation 
function and T is the time-ordering operator. 


Now we try to simplify and According to 

the Magnus expansion [^, a general time-evolution op¬ 
erator can be expanded as 


U{t) = re“*^o = exp 



(DIO) 


idip(^) 


exp 


\ 2 (w ')2 ) ■ 


2. Quantum decoherence model 


(D5) 


with the first-order and second-order Magnus terms 


Qi(t) 


—i f H{t')dt', 

Jo 

/ dt2[H{ti), H{t2)]. 
^ Jo Jo 


(Dll) 

(D12) 


In the quantum decoherence model, the NV sensor spin 
decoherence is caused by the bifurcated evolution of the 
target spin conditioned on the sensor spin state [H, [H, 
[s^ . So the decoherence of the Ah NV sensor spin under 
DD control is 


L.{t) =-Tr 


(A">)' 


U 


(-t) 


N 


(D6) 


For our specific model, the first-order Magnus term is 

S'g+e*(“=‘'-“) -hS'g-e-*(“=‘'-“) 

(D13) 


^iit) = —/ fit') 


/o 




2oji 


with 


with 


= ■■■ g-i(^^o-/3/2)(t2-ti)g-dffo-t/3/2)ti^ 

Ul^\t) = ■■■ g-i(^^o+/3/2)(t2-ti)g-dffo-/3/2)ti_ 

Eq. (|D6I) can be solved exactly for any DD control, and 
in the weak coupling regime (||/3|| ^ ||-ffo||) the sensor 
coherence dip depth can be written in a general form for 
a spin-1/2 sensor 13811. In the following part, we adapt the 
derivation in Ref. to get an expression for the NV 
sensor coherence dip depth as a function of the CPMG 
pulse number. 


F(w/,t) = w/ 


/(t')e*““*'dt' 


S'! = cos^S'f -sin^S'l', 


(D14a) 

(D14b) 

(D14c) 


Note that the term XiAl^^Sl is averaged out by the DD 
control, so it vanishes in Eq. (ID13I) . In the weak-coupling 
regime, it has been shown that ^ 2 )^) fli(l) [H, [s^. 
So we can use the first-order Magnus expansion to ap- 
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proximate and as 

« exp (-ii/oi) exp j , 


(D15a) 


Ul^\t) w exp (-ii/ot) exp 


(D15b) 


So the qubit decoherence is 


m) =2Tr 




U 


(+) 


N 


COS 




2uj' 


(D16) 


At the sensor coherence dips tdip = 7r(2(7 — l)iV/w' for the 
CPMG control, F{uj'^,t) = 2N, so the sensor coherence 
dip depth as a function of the CPMG pulse number is 


■^dipW^COS 




( A.A^,iV \ 
ywe + XiA^^^j2 j 

(D17) 


For a small CPMG pulse number, Eq. (ID17I1 from the 
quantum model agrees with Eq. (ID5II as derived in the 
semiclassical noise model. But for a large CPMG pulse 
number, the coherence dips in Eq. (jD17l) can be negative, 
while Eq. (ID5I) always predicts positive coherence dips. 
Eq. (ID 171) agrees almost exactly with the exact quantum 
results 1 ^ . 
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